home *** CD-ROM | disk | FTP | other *** search
- /* rotating beam */
- options results
-
- if ~show('ports', "AMI2D") then do
- address command 'run ami2d:bin/ami2d'
- do while ~show('ports', "AMI2D")
- end
- end
-
- address ami2d
-
- call 'ami2d:rexx/plstress'
-
- if ~show('l', "rexxmathlib.library") then do
- check = addlib('rexxmathlib.library',0,-30,0)
- end
-
- h = 0.5
- L = 10.0
- E1 = 1000
- nu1 = 0.3000
-
- reset
- iso 1 E1 nu1
- do i=0 to 20
- xx = L*i/20
- node i+1 xx 0.0
- node i+31 xx h
- end
- do i=0 to 19
- quad4 i+1 1 i+1 i+2 i+32 i+31
- end
-
- dispx 1 0
- dispy 1 0
- alias 'rotx=pbc(a1,1,0,1,10*cos(a2)-0.5*sin(a2)-10)'
- alias 'roty=pbc(a1,1,0,2,10*sin(a2)+0.5*cos(a2)-0.5)'
- rotx 51 45
- roty 51 45
- exit
-